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Abstract 

Nonrclativistic QCD provides a systematic approach for inclusive decays and productions of a 
quarkonium. By taking color-octet components into account, the approach can explain the ^/-anomaly 
at Tevatron, where the measured production rate at large transverse momentum p± is in order of 
magnitude larger than the predicted with color-singlet components only. With the approach one can 
predict that the produced J/ip and ip' at large p±_ will be transversely polarized. But the prediction 
fails in confronting with experimental measurements and this generates a puzzle. We examine the role 
of spin-flip interactions in the spin density matrix of the transition of a color-octet charm quark pair 
into J /if) and tp' . These interactions will introduce new nonperturbative parameters in the spin density 
matrix. Our result shows that the impact of the interactions is always to dilute the polarization and 
can be very significant. Taking the impact into account, predictions for the polarization are more close 
to the measured than the previous predicted. The same can also be expected for the polarization of 
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Quarkonia with their rich dynamics at different energy scales provide a special place to understand 
QCD. With its heavy mass m a heavy quark Q will move with a small velocity v inside a quarkonium. 
This results in a hierarchy of energy scales m > > mv > > mv 2 and the dynamics at different energy scales 
is different. A systematical separation of effects from the dynamics at different energy scales will provide a 
systematic understanding of quarkonium physics. Indeed, nonrelativistic QCD (NRQCD) provides such a 
systematic approach for inclusive decay and production [1] . In this approach a factorization can be made 
by expanding the small v to do the separation, a power counting in v based on the hierarchy is obtained to 
determine which nonperturbative effects, represented by NRQCD matrix elements, are relevant [H [2]. A 
comprehensive review of quarkonium physics can be found in [3]. In this approach, effects of higher- Fock 
components of a quarkonium state can be taken into account systematically. Although the probability to 
find such higher-Fock components is suppressed, but the effects of these higher-Fock states can be very 
significant. This has been shown in the explanation of the ■(//-anomaly at Tevatron, where the inclusive ip'- 
production rates with large transverse momentum p± are in order of magnitude larger than the predicted 
if one only takes the main Fock state into account. By taking color-octet QQ components into account, 
the Tevatron data[4] can be explained |i5j. This is regarded as a great triumph of the NRQCD approach. 

Despite many successes of the approach, some problems still remain unsolved. Among them a crucial 
one is that the approach fails to predict the polarization of J/ip and ip' at large p± measured at Tevatron. 
At large p±, the production of J/ip and ip' is dominant by the gluon fragmentation in which a gluon 
fragments into a color-octet cc pair in 3 Si state, the color-octet cc pair is then transmitted into J/ip and ip' . 
This color-octet mechanism successfully explains the unpolarized production rate at large p± at Tevatron. 
If the spin-symmetry for charm quarks holds, the J/ip and ip' produced through the gluon fragmentation 
is always transversely polarizedpa]. Hence, this mechanism also predicts that J/ip and ip' are transversely 
polarized at large p±. There are some effects which can dilute the polarization, like feeddown from higher 
quarkonium states, higher order corrections in a s of the gluon fragmentation [7]. Including these effects 
there is still a large fraction of J/ip or ip' with transverse polarization j8| \9\ I10j. But this prediction is in 
conflict with the measured at Tevatron|llj. This is puzzling because the same mechanism explains the 
production of unpolarized J/ip and ip' at large p± but fails to explain the polarization of the produced 
J/ip and ip' . This is the well-known puzzle. 

Giving the fact that the color-octet mechanism can explain the unpolarized production at large p±, 
for solving the puzzle one needs to carefully examine the transition of a color-octet 3 Si cc pair into ip 
which stands for J/ip or ip' , especially to examine how the polarization of the pair is transmitted. This 
is the purpose of the letter. At the leading order of v a color-octet 3 Si cc pair will be transmitted into 
a ip though the electrical dipole(£'i) interaction. This interaction conserves the spin of heavy quarks, 
hence the spin of the cc pair is totally transmitted into the spin of ip. It is also known that the magnetic 
dipole (Mi) interaction violates the spin symmetry, this will make a difference between the spin of the cc 
pair and that of ip. With the standard power counting the difference is at order of v 2 , as will be shown 
later. It is known that the suppression by v 2 for charmoina is rather weak. Hence the difference can 
be significant. The standard power counting is based on the hierarchy m » mv » mv 2 ~ Aqcd- 
For charmonia, mv is already below lGeV, and it is at the same order of Aqcd- This stimulates to 
modify the power counting for charmonia[12t [T3l HH] under the hierarchy m >> mv ~ Aqcd- in the 
modified power counting dynamical gluons are only with one typical momentum whose components are 
all at order of Aqcd ~ ttiv. With this power counting, the E\- and M\ interaction are at the same order 
Aq CD or v 2 . Since the effect of spin- flip interactions is weakly suppressed or is not suppressed, the spin 
of a color-octet 3 Si cc pair will be not transmitted into the spin of ip completely. We will examine the 
role of the spin-flip interactions with these two power-counting methods and study their impact on the 
polarization of ip' produced at Tevatron. It turns out that the impact is significant and always to dilute 
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the polarization of tp 1 . Hence the puzzle can be answered, at least at certain level. 

We consider an inclusive ip production through a color octet cc quark pair in 3 Si state. Assuming the 
produced tp is measured with the momentum p, we define the rest frame of ip by a Lorentz boost from 
its moving frame. The produced ip is polarized with the polarization vector e* in the rest frame. After 
decomposing Dirac- and color indices in NRQCD factorization, the contribution to the differential cross 
section can be generally written as: 

da[ 3 s{ 8) ] =H ij T ij (e,e*,p), (1) 

where fly is the 3x3 spin density matrix for producing a 3 Si cc pair in color-octet. This matrix can be 
calculated with perturbative QCD and contains all kinematical information. For production at a hadron 
collider, it can be written as a convolution with parton distributions and perturbative functions. The 
matrix Ty is the spin density matrix for the transition of the 3 Si cc pair in color-octet into a polarized 
ip and contains all nonperturbative information of the transition. It is defined as: 

T^(e,e*,p) = ^(0|x t TW4 c lV'(e) +*>MO +^lW t TV X |0). (2) 
x 

where the matrix element is defined in the rest frame and all NRQCD fields are at the origin of 
the space-time. The field ip(x) annihilates a c(c)-quark respectively. L is a gauge link defined as 
Pexp (— ig s J °° d\n ■ G(Xn)) with the gauge field in the adjoint representation of SU(3). n is determined 
by the moving direction of ip, it is n M = (1, — p)/\/2 with p = p/|p|. This gauge link should be added, not 
only because of the gauge invariance of Ty, but also because the completeness of NRQCD factorization 
at two-loop level, as shown by recent works [T3]. This gauge link has no effect on NRQCD factorization 
at one-loop level [l5| [IB]. 

The matrix can be decomposed as: 

Tij(e, e*,p) = Sij (e • e*a + e p e* pm) + (e^e* + e-,e*)c 

+ [(eiPj + ejPi) e* • p + (e <-> e*)] ci + pipje ■ e*c 2 H , (3) 

where • • • denote the part which is anti-symmetric in e and e* or in i and j. This anti-symmetric part is 
irrelevant here. In NRQCD the leading order interaction has the heavy quark spin symmetry. With the 
symmetry one has ao,i = c\$, = 0, i.e., only cq is not zero. This implies that ip will have the same spin as 
the color-octet cc pair does. However, the symmetry holds only approximately. Hence, in general those 
coefficients beside cq are not zero. With Eq.(3) we have the total cross-section da tot and the cross-section 
dcTi with longitudinal polarization: 

da[ Si >)tot = HijSij [2c + 3a + ai] + PipjHij [4ci + 3c 2 ] , 
da[ 3 s[ 8) ] L = HijSij [a + ai] + PiPjHij [2c + 4cy + c 2 ] . (4) 

For Tevatron, one can find that the quantity pipjHij is suppressed by relatively to the quantity 

Hijdij when p±_ is large. Therefore, for large p± the first term in the above two cross sections is dominant. 
If ao,i are zero, dai is suppressed by l/p\ in comparison with dotot and leads to the prediction that the 
produced ip will be transversely polarized at large p±. But, this is not observed in experiment [IT]. It 
should be noted that the first term can also be factorized with the gluon fragmentation function into ip 
and the function can be studied with NRQCD factorization 1 1 5 1 [T7] . 

From the above discussion, it is clear that dai will not be suppressed by l/p\, if one takes spin-flip 
interaction into account, i.e., if those coefficients beside cq are not zero. Unfortunately, these coefficients 
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are unknown yet. The corresponding coefficients of the matrix element relevant to inclusive ■(/'-decays 
has been studied with lattice QCD[18j. It has been found that the effect of spin- flip interactions is 
small. It should be noted that the matrix element for decays is different than Tj,- for productions and 
the renormalization effect is not taken into account in [18] . This effect can be difficult to study because 
one has operators which are power-divergent (see e.g., in [IH]). Although one does not know ao,i and ci^, 
one can determine their relative importance by a power counting in v, where one replaces the state ip in 
Tij with a free cc pair in a color-singlet 3 Si state and calculates with perturbative theory to extract 
those coefficients. In perturbative theory the intermediate state X in Tij should consist of two gluons 
at least. Taking the leading order interaction of NRQCD, one easily finds that cq is at order of t> 4 and 
other coefficients are zero as expected. To have nonzero ao, one of the two gluons must be emitted or 
absorbed by a heavy quark with spin-flip interaction. Taking the M\ interaction into account, one finds 
that ao is power divergent when the energy of the gluons becomes large. To have nonzero ai, c\ and C2 
an extra gluon exchanged between a heavy quark and the gauge link in Tj is required. They are also 
power divergent. It should be noted that in NRQCD the dynamical freedoms of gluons with the energy 
larger than the heavy quark mass are integrated out. Hence those gluons in the intermediate state can 
only have energies below the heavy quark mass. With a cut-off scale A s for the energy of the state X we 
have then: 

a ~ a 2 s (A s )A^, a\ ~ c\ ~ c 2 ~ a 3 (A s )A^. (5) 

The scale A s should be understood as a characteristic scale for energies of those intermediate gluons. 

If we assume that gluons with the energy at order of mv are dominant in the transition, one should 
take A s = mv. In this case the strong coupling a s should be taken as at order of t>[TJ [5J. With A s = mv 
we have 

ao ~ v 6 ai ~ c\ ~ C2 ~ v 7 . (6) 

Prom the above ao/co is suppressed by v 2 and the ratio of other coefficients is suppressed by v 3 . For 
charmonia corrections at the next-to-leading order of v are generally not small, they can be numerically 
of the same size of the leading contributions in different processes [20J. Hence those coefficients can have 
a significant impact on theoretical predictions. 

The power counting given above is the standard one. As discussed before, a modification of the power 
counting can be needed because mv is close to Aqcd for charmonia[12l [T3| fT4"] . With the modification 
a s (A s ) with A s = mv ~ Aqcd should be taken as 1. Then we have: 

£°„£l„*„^~ 0(l). (7) 
c c c c 

Therefore the coefficients ao,i and ci^ can be at the same level of importance as Co in the modified power 
counting. 

We will take the spin-flip interactions into account to predict the polarization of ip' at Tevatron. At 
Tevatron ip' can be produced through a color-singlet 3 Si cc pair and a color-octet cc pair in a 1 Sq, 3 Si 
and 3 Pj state. Beside the color-octet 3 Si state, the spin-flip interactions can also affect the transitions 
of other states. Since the color-octet 3 Si state gives the dominant contribution at large p±, we neglect 
the spin-flip interactions in contributions from other states. Summing all contributions the cross-section 
can be written as: 

da = da{ 3 s{ 1] ] + da^S^} + aV[ 3 Pj 8) ] + da[ 3 s{ 8) }. (8) 

Among the first three contributions each one is proportional to a corresponding NRQCD matrix element, 
which is (Of ( 3 Si)), (0§ ( 1 5'o)) or (Og ( 3 Pj)). These matrix elements are determined by fitting the 
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measured cross section of unpolarized ifi' at certain level. For the contribution da^s[ ], the matrix Hij 
can be extracted from the known results in We also calculated all contributions in the above and find 
an agreement with those in [9J. Without the spin-flip interactions dcTpsj 8 ^] is proportional to the matrix 
element (Og ( 3 <Si)). Taking the interactions into account, one can see that the determined (Of ( 3 Si)) is 
approximately equal to 3(2cq + 3ao + a±) from Eq.(4). 
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Figure 1: The predicted a as a function of p±. The 4 bands from the top to the bottom are corresponding 
to /o = 0, f 4 , v 2 , 1 respectively with v 2 = 0.3. The area of the bands is obtained by varying the unknown 
parameter x from to 1. In the left diagram f% is taken as v /o- In the right one /i is taken as fo 

To predict the polarization we make the ansatz that a\ = c% = C2 and introduce two parameters 
/o = oo/co and f\ = cl\/cq. These two parameters can be from to order of 1. The polarization is 
predicted with the parameter a as a function of p±, which is defined as: 

a _f da t ot s d(TL \ I ( do tot | da L \ 
V dp± dp ± J I V dp±_ dp±J ' 

If a = 1, the produced ip' is transversely polarized. If a = — 1 the produced tp' is longitudinally polarized. 
For numerical predictions we use the central values of those matrix elements used in [10] . These central 
values are determined by taking CTEQ5L parton distributions [21]. We will also use CTEQ5L parton 
distributions. With updated parton distributions these matrix elements need to be re-determined. We 
take the charm quark mass m = 1.5GeV and the energy scale of fi as fi 2 = 4m 2 + p\. It should be 
noted that for (of^So)) and (of ( 3 P )) only the combination M 3A = (of^So)) + 3A(of { 3 P )/m 2 ) 
is determined. We introduce x = (Og ( 1 »S'o))/M3.4 and x can be from to 1. For different values of /o 
and f\ we obtain the numerical results for a at Tevatron given in Fig.l. 

In Fig.l. the experimental data with errors are taken from [llj . From the figure the produced ip' will 
be dominantly with transverse polarization at large p± , if one does not take the spin- flip interaction into 
account, i.e., /q = /i = 0. Increasing /o and f\ from 0, a will be decreased. If one takes /o and fi at 
the order of 1, a is close to 0. This implies that the ip 1 is unpolarized. It should be note that /o and 
fi can be larger than 1. It is clear that the spin- flip interactions can affect the polarization significantly 
and a is more closer to the measured from experiment than that without these interactions. It should be 
emphasized that at large p± a is reduced by 50% if one takes /o = v 2 as suggested by the standard power 
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counting. If the power-counting in Eq.(7) is correct, there will be a large fraction of produced ip' with 
longitudinal polarization and the fraction can be larger than the fraction with transversal polarization. 
Although the large experimental errors prevent from an estimate of /o and fx, but one can see that 
the data favors /o ~ 1. This is unexpected with the standard power counting, but consistent with the 
power counting in Eq.(7). Although we can not determine /o theoretically, but we can examine other 
experimental data beside those from Tevatron to see if /o ~ 1 is allowed. The polarization of ip' has been 
measured from inclusive -B-decay into ip' at CLEO with a = — 0.45 ± 0.30[22j. We have made an analysis 
by adding the effect of spin- flip interactions to the relevant prediction in [23] . We find that the measured 
a at CLEO does not exclude the possibility of /o ~ 1. It should be noted that the parameter a of ip' has 
been only measured at Tevatron and CLEO. Based on our work, we can conclude that the experimental 
data do not exclude the possibility of the large spin-flip effect. 

In this letter we do not consider the case with J/ip, because the production of J/ip is more complicated 
than that of ip' . In this case the J/ip can be produced not only directly from a parton scattering, but 
also from decays of higher quarkonium states produced directly. These states can be tp' and \ c j. In 
order to understand the impact of spin-flip interactions for polarizations of J/ip, one needs to understand 
the same for those higher quarkonium states, like the case with ip' studied here. But, in general one 
can expect that the spin- flip interactions will decrease the polarization of the produced J/ip at large p±. 
The study of J/ip^s polarization, including the analysis of transitions of various cc pairs into ip with the 
spin-flip interactions and the polarization in the inclusive -B-decay will be published elsewhere. 

To summarize: We have shown that the spin-flip interactions can have a significant impact on the 
transition of a <Si cc pair into a J/ip or ip'. The impact can be parameterized by introducing new 
parameters in the transition matrix T. If the heavy quark spin symmetry holds, the matrix has only 
one parameter, denoted as cq. Because the charm quark is not heavy enough, or v is not small enough, 
these new parameters are not small in comparison with cq and can be at the same size of cq. These 
parameters can significantly reduce the polarization parameter a of ip' , measured at Tevatron. Including 
the effect of the spin-flip interactions, the predicted polarization of ip 1 will be close to the measured and 
the possibility of the large spin-flip effect is not excluded by known experiments. This provides a way to 
solve the puzzle discussed at the beginning. It is also expected that the spin-flip interactions will change 
the polarization J/ip significantly and can provide a solution to solve the puzzle of J/ip. 
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